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A novel singular gauge transformation is developed for quasiparticles in the mixed state of a strongly 
type-II superconductor which permits a full solution of the problem at low and intermediate fields, 
H c \ < B <C H C 2- For a periodic vortex lattice the natural low-energy quasiparticle states are 
Bloch waves rather than Landau levels discussed in the recent literature. The new representation 
elucidates the physics considerably and provides fresh insights into the spectral properties of such 
systems. 



In conventional (s-wave) superconductors the single 
particle excitation spectrum is gapped and, consequently, 
no quasiparticle states are populated at low tempera- 
tures. The situation is dramatically different in uncon- 
ventional superconductors which exhibit nodes in the 
gap. These lead to finite density of fermionic excitations 
at low energies which then dominate the low-temperature 
physics. Among the known (or suspected) unconven- 
tional superconductors are high-T c copper oxides, organic 
and heavy fermion superconductors, and the recently dis- 
covered Sr2Ru04. Understanding the physics of the low- 
energy quasiparticles in the mixed state of these uncon- 
ventional superconductors is an unsolved problem of con- 
siderable complexity. This complexity stems from the 
fact that (i) in the mixed state superconductivity coex- 
ists with the magnetic field B and the quasiparticles feel 
the combined effects of B and the spatially varying field 
of chiral supercurrents, and (ii) being composite objects, 
part electron and part hole, quasiparticles do not carry 
a definite charge. The corresponding low-energy theory 
therefore poses an entirely new intellectual challenge |Q , 
which is simultaneously of considerable practical interest. 

The initial theoretical investigations were based on nu- 
merical computations semiclassical approximations 
H and general scaling arguments [Q. More recently 
Gorkov and Schrieffer |^] made a remarkable prediction 
that in a d x 2_ y 2 superconductor at intermediate fields 
H c i <C B <C H C 2 the quasiparticles will form Landau 
levels (LL) with a discrete energy spectrum 

E n = ±htu H \/n, n = 0, 1,..., (1) 

where u>h = 2\Jlu c /S.q/Ti, with uj c = eB/mc being the cy- 
clotron frequency and Ao the maximum superconducting 
gap. Based on a somewhat different reasoning Anderson 
U later arrived at a similar conclusion and argued that 
LL quantization could explain the anomalous magneto- 
transport in cuprates. Janko j?J proposed a direct test of 
Eq. (|l|) using the scanning tunneling spectroscopy. 

The concept of LL quantization here is quite differ- 
ent from the one in conventional superconductors, where 
nodes in the gap arise as a result of the center of mass 
motion of pairs in strong magnetic fields near H C 2 
The physics of Eq. fll]) is based on the picture of a low- 



energy quasiparticle Larmor-precessing in weak external 
magnetic field along an elliptic orbit of constant energy 
centered at the Dirac gap node in the fc-space . This 
motion corresponds to a closed elliptic orbit in the real 
space and the quantization condition (|l|) follows from de- 
manding that such an orbit contains n quanta of mag- 
netic flux. This argumentation neglects the effect of spa- 
tially varying supercurrents in the vortex array, which 
were, however, recently shown by Melnikov Q to strongly 
mix the individual Landau levels. 

In this paper we formulate a new approach to the prob- 
lem which treats the effects of the magnetic field and 
supercurrents on equal footing. As a result the physics 
becomes transparent and for periodic vortex arrays in 
the intermediate field regime the low-energy theory can 
be solved in its entirety. Our principal result is that the 
natural low-energy quasiparticle states are Bloch waves 
of massless Dirac fermions and not the Landau levels dis- 
cussed above. 

At the heart of our approach is the observation that 
the collective response of the condensate to the external 
magnetic field on average exactly compensates its effect 
on the normal quasiparticles. More formally, the phase of 
the superconducting order parameter, A(r) = Aoe l ^ r \ 
acts as an additional "gauge field" coupled to the quasi- 
particles. In the vortex state 4>{r) is not a pure gauge: 
V x V0(r) = 2nz J2i <5( r — Hi) where {Ri} denotes vor- 
tex positions. From the vantage point of a quasiparticle 
the singularities in V x V0 act as magnetic half-fluxes 
concentrated in the vortex cores with polarity opposing 
the external field. Flux quantization ensures that on av- 
erage this "spiked" field exactly cancels out the external 
applied field B. In the mixed state, the quasiparticle 
therefore can be thought of as moving in an effective 
field B ff which is zero on average but derives from a 
vector potential that is highly nontrivial. The nature 
of the phenomenon is purely quantum-mechanical, and 
is closely related to the Aharonov-Bohm effect: classi- 
cal charged particle would be completely unaffected by 
the spiked field because the singularities occupy a set of 
measure zero in the real space. 

Our solution consists in finding a gauge in which the 
Hamiltonian manifestly displays the physical property 
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described above. In such a gauge the fermionic excita- 
tion spectrum can be found by band structure techniques 
suitably adjusted to the "off-diagonal" structure of the 
theory. Besides revealing the nature of the low-energy 
quasiparticles, this representation leads to new insights 
into the physics of the mixed state. One surprising find- 
ing is that in a perfectly periodic vortex lattice the orig- 
inal Dirac nodes survive the perturbing effect of a weak 
magnetic field. 

We now supply the details. Quasiparticle wavefunc- 
tion ^/ T (r) = [u(r),i;(r)] is subject to the Bogoliubov-de 
Gennes equation H$ = E^>, where 



n 



He A 

A* -Ht 



(2) 



with He — 2m(P — f A) 2 — ep and A the d-wave pair- 
ing operator. Following Simon and Lee Q we choose the 
coordinate axes in the direction of gap nodes, in which 
case A = Pp 2 {p x , {py, A(r)}}, where pf is the Fermi mo- 
mentum, p = —iKV, and curly brackets represent sym- 
metrization, {a, b} = ^(ab + ba). In high-T c cuprates 
it is natural to concentrate on the low to intermediate 
field regime H c \ < B <C H c2 , where the vortex spacing is 
large and we may assume the gap amplitude to be con- 
stant everywhere in space, A(r) ~ Aoe l< ^ r - ) . Under such 
conditions the magnetic field distribution is described by 
a simple London model and superfluid velocity, defined 
as v s (r) = ^;(§V0 — |A),_ can be written in terms of 



vortex positions {R;} as [|10[ 



v s r = — 

m 



d 2 k ik x z ^ 
(2^)2 A- 2 + P ^ ' 



ik-(r-Rj) 



(3) 



where A is the London penetration depth. 

In order to diagonalize (||) it is desirable to remove the 
phase factors e^ r ' from the off-diagonal components of 
H. This is accomplished by a unitary transformation 



H -> U^HU, U 



g*0e(r) 







,-ic/> h (r) 



(4) 



where </> e (r) and </>/t(r) are arbitrary functions satisfying 

e (r)+^(r)=0(r). (5) 

Eq. ([|) can be thought of as a singular gauge transfor- 
mation since it changes the effective magnetic field seen 
by electrons and holes. We now discuss three specific 
choices for the functions <p e and <ph- 

The most natural choice satisfying (|^) is the symmetric 
one, namely </> e ( r ) = ^Ti(r) = </>(r)/2, resulting in Q 



H s = 



J—(v + mv 

A 



l PxPy 



PxPy 



J-(p-mv s ) 2 + e F 



This particular gauge makes the Hamiltonian very sim- 
ple but unfortunately is not very useful because, as noted 



by Anderson |6) and in a different context by Balents et 
al. |]l2|| , the corresponding transformation (Q) is not sin- 
gle valued. To see this, consider the situation on en- 
circling the core of a vortex: (j) winds by 2-7T but <p e 
and 4>h pick only a phase of 7r, causing U to have two 
branches. Consequently, one is forced to diagonalize Hs 
under the constraint that the original wavefunctions are 
single valued. Clearly, this is a difficult task. Never- 
theless, the symmetric gauge reveals the physical essence 
of the problem: formally, v s enters Hs as an effective 
vector potential, which corresponds to an effective mag- 
netic field B c ff = - j(V x v s ) / B. It is easy to show 
from Eq. (||) that B c ff vanishes on average |l3|], i.e. that 
(V x v,) = 0, where angular brackets denote the spatial 
average. Aside from the single-valuedness problem, the 
low-energy physics described by Hs is that of a quasipar- 
ticle in zero average magnetic field. The external field is 
compensated by the array of magnetic half-fluxes giving 
rise to a non-trivial vector potential with the periodicity 
of the vortex lattice. 

To avoid the problem of multiple valuedness Anderson 
||, suggested taking (j> e {r) = ^(r) and 0h( r ) = 0. This 
leads to a Hamiltonian of the form [pTl 



Ha 



cr-(p+ - A + 2mv s 
D 



D 

5±-(p + ^A) 2 



with D = j^-(p x + ^A x +mv sx )(p y + ^A y +mv S y). In this 

representation one could consider neglecting in the first 
approximation the v s terms, on the grounds that they 
represent a perturbative correction ||. At low energies, 
expanding all the terms in Ha to leading order near the 
nodes, the Hamiltonian becomes that of massless Dirac 
fermions in a uniform magnetic field B = V x A with 
the energy spectrum || given by Eq. ([!]). The effects of 
supercurrent field v s can in principle be treated pertur- 
batively. This is, however, difficult in practice, because of 
the massive degeneracy of the Landau levels and also be- 
cause v s (r) is not a small perturbation (it diverges as ljr 
at the vortex cores) and will lead to strong LL mixing M . 
In the absence of a reliable scheme to incorporate v s (r) , 
the physical picture of Larmor precessing quasiparticle 
that leads to Eq. (|l]) appears incomplete. 

We now introduce a new singular gauge transformation 
that combines the desirable features of the two transfor- 
mations discussed above but has none of their drawbacks. 
Consider dividing vortices into two distinct subsets A and 
B, each containing an equal number of vortices. Now de- 
note by 4>a{y) the phase field associated with vortices in 
the subset A, with the analogous definition of 4*b (r) . The 
choice e (r) = </>a(i")> 'Ph (r) = </>s(r) clearly satisfies the 
condition dq) and the corresponding transformation U is 
single valued. The resulting Hamiltonian is 



H 



N 



i-(p + mvf) 2 - £f 



D 



D 



CF 



with D = $\p x + f{v? x - vf x )][Pv + f K ~ vf y )] and 
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FIG. 1. (a) Two sublattices A and B of the square vortex 
lattice and the unit cell used in the numerical diagonalization 
of Eq.(H). (b) The corresponding magnetic Brillouin zone. 



^ = -(KV<f> M --A), (m = A,B. 
m c 



(6) 



As long as the physical field B remains approximately 
uniform, v^(r) satisfy equations similar to (||) with Ri 
replaced by R^ and an overall prefactor 2. 

It is easy to verify that v% correspond to zero effective 
field, i.e. that (V x v^ 1 ) = 0. Evidently, there is no reason 
to expect LL quantization in the system. This property 
becomes more transparent if we focus on the low-energy 
excitations. By linearizing TL^ in the vicinity of the four 
nodes as described in Ref. j| we obtain H N ~H +H', 
where 



Ho 



VFPx VAPy 
VAPy -VFPx 



is the free Dirac Hamiltonian and 

\v A {v 



\va{v 
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(7) 



(8) 



is the vector potential term, vp is the Fermi velocity and 
va = Ao/pf denotes the slope of the gap at the node. 

Our considerations so far have been completely general 
and apply to arbitrary distribution of vortices. In the fol- 
lowing we illustrate the utility of the new Hamiltonian by 
finding the excitation spectrum in a periodic square vor- 
tex lattice. With minor modifications the same approach 
can be generalized to an arbitrary periodic lattice, such 
as e.g. triangular. We take A and B subsets to coincide 
with the two sublattices of the square vortex lattice as 
illustrated in Figure 0(a) . Expanding the quasiparticle 
wavefunction in the plane wave basis \I/(r) = 5' q e lq ' r 
we arrive at an equation of the form 



Wo(q)*q + ^W'(K)* q+K = 



(9) 



K 



Because of its periodicity, Ti' only has non-vanishing 
Fourier components at the reciprocal lattice vectors 
K = =j-(m x ,m y ), where (m x ,m y ) are integers and d = 

y2$o/B is the size of the unit cell. Aside from the 2x2 
matrix structure, Eq. (Q) is a standard Bloch equation 
which we solve by numerical diagonalization for q vec- 
tors in the first MBZ, sketched in Figure |l|(b). 

As long as A ^> d the results are independent of A and 
the only free parameter in the low-energy theory is the 




FIG. 2. Upper panel: lowest energy bands near a single 
node for isotropic case an = vf/va = 1- Dashed lines repre- 
sent the unperturbed Dirac spectrum of Ho while solid lines 
reflect the effect of the periodic potential Ti! , Lower panel: 
the corresponding DOS. Dashed line represents the thermally 
broadened DOS resulting from LL spectrum of Eq. (N). 



Dirac cone anisotropy old = vf/va- The band structure 
for the isotropic case, «j) = 1 is presented in Figure |^. 
When compared to the unperturbed Dirac bands the pe- 
riodic potential has precisely the expected effect of open- 
ing up band gaps at the MBZ boundaries. The surpris- 
ing finding is that the magnetic field does not destroy the 
original nodal point, but merely renormalizes the slope 
of the dispersion. This finding can be understood as a 
consequence of the exact electron-hole symmetry of the 
linearized Hamiltonian The associated DOS van- 

ishes at the Fermi level, in contrast to thepeak expected 
from the LL scenario (cf. Figure g). The peaks 

which appear in DOS are van Hove singularities related 
to the band structure and have nothing to do with LL 
spectrum of Eq. ([!]). Observation of these peaks is a 
challenge to the experimental community. 

In Figure |we display the band structure for old = 20, 
a value perhaps more relevant for the optimally doped 
YBa2Cu307_5. The striking new feature is the forma- 
tion of additional lines of nodes on the Fermi surface [see 
also Figure [fjb)] which give rise to a finite DOS at the 
Fermi surface. This structure can be understood by con- 
sidering the old — > oo limit Q . The lines of nodes first 
appear at old ~ 15 and with increasing nodal anisotropy 
their number increases. This finding is consistent with 
the prediction of a finite residual DOS based on the semi- 
classical approach M , which is expected to be valid when 
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FIG. 3. The same as Figure | with a D = 20. 

Ao <C e_F, or equivalcntly up 3> 1. Qualitatively similar 
results are found for different orientations of the vortex 
unit cell with respect to underlying ionic lattice. 

In conclusion, we have shown how the collective su- 
perfluid response of a superconductor ensures that the 
effective magnetic field B e fj seen by a fcrmionic quasi- 
particle (and distinct from the physical field B) is zero 
on average, even in the vortex state. The physics of a 
low energy quasiparticle in a superconductor with gap 
nodes is that of a massless Dirac fermion moving in a 
vector potential associated with physical supercurrents 
but zero average magnetic field. For a periodic vortex 
lattice the appropriate description is in terms of familiar 
Bloch waves. A mathematically equivalent description 
could be given (in a different gauge) in terms of Lan- 
dau levels strongly scattered by supercurrents ||. 
The fact that no trace of LL structure remains in the ex- 
act spectrum of excitations suggests that the former is a 
more useful starting point. The LL quantization remains 
a domain of relatively high fields B. Our conclusions 
are corroborated by the absence of LL spectra in the nu- 
merical computations Q as well as in the experimental 
tunneling data on cuprates and are consistent with 
scaling arguments given previously ||]4] . 

In the present study we have focused on the leading 
low-energy, long-wavelength behavior of the quasiparti- 
cles as embodied by the linearized Dirac Hamiltonian (0- 
|8|). In real materials and at higher energies our results 
may be modified by the corrections to the linearization, 
electron-hole asymmetry possible internode scattering, 
ionic lattice effects and the vortex core physics. These 
issues, as well as a more rigorous discussion of the singu- 



lar gauge transformation, are best addressed within the 
framework of a tight binding calculation to be reported 
in a forthcoming publication. 

We emphasize that the central idea of this paper, i.e. 
that upon proper inclusion of the condensate screening 
the quasiparticles experience effective zero average mag- 
netic field, is completely general and robust against any 
effects of short length scale physics. Consequently, our 
method is applicable to any pairing symmetry and arbi- 
trary distribution of vortices and could be useful for un- 
derstanding the physics of vortex glass and liquid phases, 
as well as the zero-field quantum phase-disordered states 
such as the nodal liquid [ ^2[ . We expect that disorder 
in the vortex positions will smear the structure apparent 
in Figures land | resulting in smooth DOS. Of obvious 
interest are implications for the quasiparticle thermody- 
namics, transport and localization properties in statically 
disordered or fluctuating vortex arrays. 
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